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1. Introduction 

Self-avoiding walk models are of interest both as models with which to understand the 
thermodynamic behaviour of polymers in solution, and as models of theoretical interest 
in thier own right pi El IHj . 

The thermodynamic behaviour of a linear polymer in dilute solution is determined 
by the competition between an attractive interaction, due to the difference between 
the monomer-monomer affinity and the monomer- solvent affinity, and an effective 
repulsive interaction (the excluded-volume interaction) which results from the loss in 
entropy caused by bringing together two segments of polymer. At high temperature 
the excluded-volume interaction dominates, leading to the good solvent phase, where 
the polymer adopts an open configuration, well interpenetrated by the solvent. At low 
?temperatures the monomer-monomer affinity dominates and the polymer collapses into 
a dense ball, and in practice precipitates from solution. This is the bad solvent phase. 
These two phases are separated by the 9-point transition jH Ej. 

The thermodynamics of a polymer in solution is well modelled on the lattice by 
the so called G-point model|ni Ej, which consists of a self-avoiding walk on a regular 
lattice with the inclusion of attractive interactions between non-consecutively visited 
nearest-neighbour lattice sites. This model displays behaviour in excellent agreement 
with experimental results in both two and three dimensions. The critical behaviour of 
the "on-lattice" model is in agreement with the results of equivalent "off-lattice" models, 
indicating that the lattice does not affect the critical behaviour. 

The interesting feature of the B-point model is that it describes the thermodynamic 
behaviour of a large class of linear polymers, independently of the detailed chemistry 
of the polymer chain. This may be understood as follows: the real polymer chain, 
if long enough, behaves as a chain of spherical globules, of comparable in size to the 
persistence length. The detailed chemistry simply defines the size of the globules, and 
the large scale behaviour of the chain is described by a random walk, constrained by the 
excluded volume condition (globules repel each other due to entropic repulsion). The 
globules interact via effective attractive interactions, incorporating a variety of different 
microscopic interactions. The monomers in the model system then corresponds to a 
large number of real monomers on the scale of the real polymer chain. Bearing this 
in mind, when going over to a lattice model, the decision to think of the 'monomers' 
sitting on the lattice sites or the lattice bonds is arbitrary. The choice of placing the 
interactions between the the bonds would be, in principle, just as reasonable. 

In this article we present a transfer matrix study of the bond-interacting self- 
avoiding walk model, and show that the phase diagram differs in many important 
ways from the standard 0-point model, displaying richer behaviour. This model has 
been previously studied using extended mean-field type calculations |H1 El HD] • Different 
methods used by the different authors resulted in two different proposed phase diagrams. 
Support for the conjecture due to Stilck and co-workers [HI E] (summarised in figure [l]) 
is provided by Machado et al [llj, who studied the model also using transfer-matrix 




Figure 1. Phase diagram proposed by Stilck and co-workers [HI ITT] . The lower Hne 
is a critical line in the self-avoiding walk universality class, conjectured to terminate 
in a critical end point, CEP. The upper line is conjectured to made up of a critical 
transition line and a first order transition line separated by a tricritical point, TCP. 
Phase I is the finite walk phase, phase II an isotropic dense phase, and phase III an 
anisotropic crystalline phase. 



methods. Our results differ in important details, and lend support to the phase diagram 
proposed by Buzano and PrettiJU]- This phase diagram is shown schematically in 
figure 13 We shall attempt to explain the reasons for these discrepancies. 



2. The model and the transfer matrix method 



The model studied in this article is the self-avoiding walk on the square lattice. A 
chemical potential /x, or equivalently a fugacity k = exp(— is associated with each 
step of the walk. An additional attractive energy —e is introduced every time two steps 
of the walk are parallel to each other across the face of a lattice plaquette, see figure El 
The thermodynamics of this model may then be studied in the grand-canonical ensemble 
by introducing the partition function: 

Z=J2 ^''exp{Njf3e), (1) 

walks 

where N is the number of steps in the walk, Nj is the number of interactions, and 
f3 = 1/kT as usual. It is convenient in what follows to set £ = 1, which is simply 
equivalent to a change in the temperature scale. 

Once the partition function has been found, many thermodynamic quantities follow. 
The (dimensionless) free energy per site is defined as: 




Figure 2. A conjectured schematic phase diagram, summarising the results found in 

this article. The phase diagram splits into three phases: (I) a finite walk phase {p = 0), 
(II) an isotropic dense walk phase (liquid phase) and (III) an anisotropic dense walk 
phase (crystalline phase). The transition from phase (I) to phase (II) is critical in the 
self-avoiding walk universality class for (3 < j3i and first-order between [3i and (32 (in 
analogy with the standard Q point model) . The transition from phase (I) to phase (III) 
is first order. The transition occurring at j3i is a tricritical point which we conjecture 
to be in the © universality class, whilst the transition at /?2 is a critical end point. The 
transition from phase (II) to phase (III) is conjectured to be in the BKT universality 
class. 



V V 




Figure 3. Example of a walk configuration, showing example column states. The 
nearest- neighbour interactions are indicated using dashed lines, whilst the thin lines 
show the connectivities that need to be taken into account in the column state. 
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where Q denotes the number of sites on the lattice. The density, p, of the walk on the 
lattice and the energy per lattice site follow: 

(N) Of 

Further derivatives give the response functions, the susceptibility and the specific heat. 

In general there is no exact analytical result for the partition function, and so the 
thermodynamic quantities must be calculated either by simulation (for example using 
the Monte-Carlo method) or some approximation scheme must be found enabling the 
determination of the partition function. In this article we shall use the transfer matrix 
method, where the exact partition function can be calculated for a lattice of finite width 
L and of infinite length O dlllllS US El . The width of the lattice can then be 
varied, and the results extrapolated to infinite lattice width. 

The standard way of considering the problem is as follows: define the restricted 
partition function Zx{Co, C^) as the partition function for a portion of the lattice between 
a; = and x. The walk has a column state Cq at the origin and Cx in column x. One 
may then write the following recursion relation: 

Zx+liCo, Cx+l) = ^ Zx{Cq, Cx) T{Cx, Cx+i), (5) 

where T{Cx, Cx+i) is the additional Boltzmann weight to add column x + 1 in 
configuration Cx+i next to column x in configuration Cx- This forms a (transfer) matrix. 

That this recursion should be valid for a spin system is fairly straightforward, since 
the interactions are all local. For a polymer model it is less clear that this should be 
possible, since one has to take into account non-local factors, most notably to ensure 
that the partition function describes only one chain, without the formation of "orphan" 
loops. This is done by appropriately defining the column states Cx- For the self- 
avoiding walk problem with no added interactions (i.e. e = 0), it is sufficient to define a 
column configuration by the arrangement of horizontal bonds in the column along with 
information about the connectivities between the bondsfSJ El E]i i-e- information 
about which pairs of horizontal bonds are connected by polymer loops to the left (taking 
X as increasing towards the right). See figure El for an example of how to define the 
column states for the bond- interacting model. By successive application of (0), one 
finds 

ZxiCo,Cx) = {Co\T-\Cx)- (6) 

Since we are interested in equilibrium thermodynamic behaviour, and we shall be taking 
the limit x — oo, the choice of boundary conditions is arbitrary. Choosing periodic 
boundary conditions in the x direction results in Zx = Tr T^, giving: 

= E (7) 
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terms of the eigenvalues Aj of the transfer matrix T. In general the largest eigenvalue is 
non-degenerate, and the sum is dominated by this largest eigenvalue, Aq, giving, in the 
limit X — oo, 

/ = ^logAo, (8) 

where L is the width of the lattice strip. 

Having calculated the free energy, it is now possible to calculate the density, the 
energy and other local quantities by differentiation. However, since the eigenvalues of 
interest will be determined numerically, it is better to avoid differentiation whenever 
possible; the eigenvalues themselves may be calculated to arbitrary precision, but the 
numerical differentiation magnifies round-off errors, leading to a substantial loss of 
numerical precision. This is not a problem, however, since in the transfer matrix 
framework it is easy to calculate averages over local quantities directly[l8]. To see 
this, it is necessary first to calculate the probability of having a given configuration Cx 
in column x. This probability is simply the ratio of the partition function restricted 
to configuration and the unrestricted partition function, which in terms of transfer 
matrices may be written: 

where M is the length of the lattice strip. 

Writing | Cx) in terms of the eigenvectors, \i) of T gives: 

E^Xf'{^\C ){C\^) 



p(g) = ,lim, ^---^-'-;^-'-^ (10) 



p{C) = {0\C)', (11) 

where the eigenvectors are normalised. The subscript x may be omitted by invoking 
translation invariance. 

The density, for example, is then found using 

sr^N{C) 
P = 2^—i—p{G) 

c ^ 

c ^ 

where N{C) is the number of occupied lattice bonds in configuration C. The 
susceptibility can then be calculated either by taking a derivative of the density, or by 
calculating directly (A^^) for the column, and hence the fiuctuation. The two methods 
give slightly different results for a finite width strip, but agree in the thermodynamic 
limit. In the present article we choose to calculate the fiuctuation directly. 

It is also possible to calculate two-point correlation functions, by considering the 
joint probability p{Cx, Cy) of having a configuration Cx in column x and Cy in column 
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y. The calculation proceeds exactly as above: 

P^^- = i^-.^ Tr{T-} ^''^ 

lim -^{Cy\T''^^-'^\C.){C.\Ty-^\Cy) (14) 



^0 \ 



(15) 



where a and b are constants and Ai is the second largest eigenvalue. This leads directly 
to the identification of the correlation length as|19j 

The correlation length diverges when the two largest eigenvalues become degenerate, 
signalling the onset of a phase transition. For many spin models there is no phase 
transition below two dimensions, and so the correlation length should not diverge 
for finite lattice widths. In the transfer matrix formulation this is seen explicitly by 
considering the Frobenius-Perron theorem, which states that for a finite-dimensioned 
matrix with strictly positive entries, the largest eigenvalue is non-degenerate PU] . In our 
model, however, there are many zero entries, and the Frobenious-Peron theorem does 
not apply. The potential degeneracy of the eigenvalues may be used to gain an insight 
into the phase diagram, and we shall return to this fact below. 

The correlation length can be used as the basis of a phenomenological 
renormalisation group analysis |2T]. The idea is simple: far from the critical point the 
correlation length is smaller than the width of the lattice, but as the critical point is 
approached the divergence of ^ is limited by the smallest dimension, L. Since we expect 
to have scale invariance when both ^ and L ^ oo, all thermodynamic quantities should 
(for sufficiently large lattice sizes) be only by a function oi ^/L. The critical point may 
then be identified with solutions of the equation: 

^ = ^. (17) 

Two lattice sizes are required to determine the estimate, reducing the number of available 
estimates for extrapolation. There is no guarantee that there is a solution to this 
equation at finite lattice widths, and it is sometimes necessary to show that the gap 
goes to zero in the thermodynamic limit. There is also no guarantee that a solution to 
this equation corresponds to a transition point. Assuming the solution is an estimate 
for a phase transition, then using the asymptotic behaviour ^ ~ \k — Kc\~'^, we may infer 
estimates for z/ at the same time: 

1 log(1^/y 
UL,L' log (L/U) 

Other quantities may be used in a phenomenological renormalisation scheme. For 
self-avoiding walk type models the density is of particular interest. Close to a critical 



-1. fl^ 
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point the singular part of the free energy is expected to scale as: 

fsiK,L) = L-''fi\K-K,\L'/n, (19) 

where d is the spatial dimension, here d = 2. Taking the first derivative with respect to 
K gives: 

Ps = L'/''-'p{\^^-K,\L'/n- (20) 

We have denoted the density as ps since this is the density related to the singular part 
of the free energy, which vanishes in the thermodynamic limit. In general there will be a 
non-zero contribution from the non-singular part of the free energy, giving the limiting 
density. In walk models, the Iow-k phase corresponds to finite length walks, and hence 
the density in the thermodynamic limit is zero. The scaling law then enables us to 
define a renormalisation scheme [22| I28j. We start by defining the scaling function: 



log {p,{k,L)/ps{k,L')) 

= ki(I7I7)) • pi) 

Since we wish to look for a family of estimators for the critical point, and taking into 
account that there are potentially severe parity effects, it is convenient to set L' = L — 2. 
Estimates of the critical point are then found by looking for solutions of the equation: 

<fL,L-2{fil) = VL-2,L-4{i^I)- (22) 

If such a solution exists, then = limi^oo'^L and z/ = limi^oo(l/ {2 + ^pL,L-2{'^L))■ 
0theT quantities will be used later, but the idea remains the same. 

Of course estimates of transition points can also be found by looking at appropriate 
susceptibilities, and identifying the transition with the peak of the function. 

Each method has its advantages and disadvantages. The phenomenological 
renormalisation group method proposed by Nightingale is usually the most asymptotic, 
but we shall argue here that some of the lines found using this method do not correspond 
to transition lines, most notably that found by Machado et al ^l] . Caution is therefore 
important, and it is wise to crosscheck results using different methods. 



3. The K—[3 Phase Diagram 

In this section, the phase diagram is mapped out in the n-[3 plane. 

It is expected that the model will have a behaviour similar to the G-point model 
as K is increased from zero: at low enough k the length of the walk is finite, and 
since the walk sits on an infinite lattice, the density is zero. The average length of 
the walk diverges along the line k = k,i{(3) , either smoothly, for small enough (3, or 
discontinuously, for large enough (3. In the first case we have a critical phase transition 
in the same universality class as the self-avoiding walk model {(3 = 0) and in the second 
we have a first order transition. In the G point model the point separating the two 
regimes is a tricritical point, with a correlation length exponent u = 4/7, different 
from the self-avoiding walk value u = 3/4. Similarly, as f3 is increased along the line 
K = Ki{(3), the interactions of the bond-interacting self-avoiding walk will drive a collapse 
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transition. In the B-point model there are no other phase transitions, but for the bond- 
interacting model extended mean-field calculations predict a phase transition between 
two dense walk phases |H1 El CD], an isotropic phase (the liquid phase II in figure |21) and an 
anisotropic phase (the crystalline phase III in figure |21). We will identify this transition 
line, K = K,2{P), with a crystallisation transition. This line will terminate on the line 
K = Ki{i3) at a value of (3 = (32- An important question to be answered is whether 
Pi = P2, i-e., does the walk, as the line k = Ki(/3) is followed, collapse then crystallise, 
or rather collapse directly to a crystalline, anisotropic, state? The results we present 
bellow tend to support the first scenario. 

As mentioned in the previous section, the correlation length diverges when the two 
largest eigenvalues become degenerate. As a result of the connectivity constraints, this 
may occur in the walk models at finite lattice width. In the Iow-k phase the density 
is zero, and the corresponding largest eigenvalue is equal to one for all values of k. 
The second largest eigenvalue steadily increases as k increases, and for some value of 
K = Kl{i3) the two eigenvalues become degenerate. For n > the walk will fill the 
lattice with a finite density. The largest eigenvalue is now bigger than one, and clearly 
does not correspond to the vacuum state. An estimate for the location of the transition, 
n = is then determined by equating the two largest eigenvalues Aq = Ai = 1. 

It turns out that the transfer matrix block-diagonalises naturally into three sectors: 
the vacuum state forms a block on its own, an odd sector (corresponding to walks with 
an odd number of links per column) and an even sector (where the walks have an even 
number of links per column). Since the vacuum state is trivial, and the eigenvalue is 
constant and equal to one, no calculation is required. Usually the calculation of the 
transfer matrix is limited to the odd sector and then the transition point is determined 
by simply setting Ai = l^Hl- Since the smallest possible increase of density on a finite 
width lattice is Ap = 1/L, this transition, which exists at finite lattice widths, will 
appear as a first-order transition, becoming critical in the infinite lattice width limit. 

This method is easy to implement, and often a good choice to find estimates of 
the critical point, since we find one estimate per lattice width L. Since the number 
of lattice widths one can calculate is limited by the rapid increase in the size of the 
transfer matrix, it is important to use estimators that require the smallest number of 
lattice sizes if we are to apply a reasonable extrapolation scheme. The disadvantage of 
this method is that the estimates of the critical point may start reasonably far from the 
true critical point. 

In what follows, we define Ai as the largest eigenvalue of the odd sector and A2 as 
the largest eigenvalue of the even sector. 

The even sector of the matrix corresponds to walk configurations that cross the 
lattice an even number of times. The whole analysis may be repeated with the even 
sector, and it would be expected that the same final results would be found, the difference 
between the walk crossing once or twice simply corresponding to a change in boundary 
conditions. For the standard 6-point model, Ai > A2 everywhere, and so the lines found 
are simply less asymptotic estimates. This is no longer true for the bond-interacting 
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Figure 4. Estimates of the zero-density/finite-density phase transition hne determined 
by setting A = 1 for the even sector (open symbols) and odd sector (closed symbols). 
The solid lines indicate the transition line estimates. The collapse transition is 
identified with the point Ai = A2 = f . 

self-avoiding walk for even lattice sizes, where for high enough /3, A2 reaches one first 
as K increases. This is shown in figured There is a point, then, on the k = Hi{f3) 
line in which we simultaneously have Ai = A2 = 1. We identify this point with the 
crystallisation transition. This is a reasonable identification if one considers that in the 
crystalline phase the walk is trying to align its bonds with a lattice direction, and so 
must minimise the number of corners. For even lattice sizes the easiest way of achieving 
this is for the walk to fold back onto itself an even number of times. The corresponding 
configurations sit in the even sector of the transfer matrix, leading to the observed 
crossover in the eigenvalues. This is of course not the case for odd lattice widths, where 
the same argument leads to an odd number of folds, with corresponding configurations 
which sit in the odd sector leading to no such degeneracy in eigenvalues, i.e. Ai > A2 
everywhere. 

By extension, the condition Ai = A2 defines a line in the k — (3 plane. When 
Ai = A2 < 1, then the largest eigenvalue corresponds to the vacuum state Aq = 1. 
Crossings of subdominant eigenvalues indicate a change of local order, and define so- 
called disorder lines |24[ |23 12^1 12Z1- However, when Ai = A2 > 1 we have a crossing of 
the largest and second largest eigenvalues, and thus have a divergent correlation length. 
This condition defines a transition line, possibly critical, in the dense walk phase. We 
identify this line with the crystallisation transition line k = K2{(3). The full phase 
diagram found using these arguments is shown in figure 

Machado et al performed a transfer matrix calculation for this model, and 
found evidence of a transition line, though at a much lower value of (3 than the transition 
line predicted here. In order to resolve this discrepancy, we look for solutions to the 
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Figure 5. Complete phase diagram found using eigenvalue crossings. The lower line 
is found using the condition max(Ai, A2) = 1 wliilst the upper line is found using the 
condition Ai = A2. 

phenomenological renormalisation group equation (jl7|) . Results for even lattice sizes 
are shown in figure IHl In determining the correlation lengths, best results are found if 
the two largest eigenvalues are used. To be sure of having the two largest eigenvalues, it 
is necessary to calculate in both the odd and the even sectors. It turns out that in the 
region of interest the two largest eigenvalues appear in different sectors, which simplifies 
their calculation. These results are complementary to those given by Machado et al 
[TT] , who limited their study to the odd sector of the transfer matrix, and mostly to odd 
lattice sizes. 

In figure ini we see the appearance of an additional line, which is consistent with the 
transition line proposed by Machado et al ^T]. 

The question now arises: which lines correspond to true phase transitions? It 
is natural to look for indicators such as response functions (fiuctuations in the order 
parameter, specific heats, etc.), which may be expected to show a singularity at the 
phase transition. A first step is to define a suitable order parameter; a good choice is 
the difference in densities of horizontal and vertical bonds, 5p = \ph — Pv\- The energy 
is maximised when all the monomers are parallel and pointing in the same direction. 
The effect of increasing temperature is to disorder this ground-state, mixing up the two 
phases, and restoring isotropy. Figure [7| shows both 6p, and corresponding fluctuations, 
for K, = 0.5, i.e. within the dense walk region of the phase diagram. The response 
function is xsp = L{{6p'^) — {Sp)^). Figure |H1 on the other hand, shows the density 
and its corresponding fluctuations. The density of interactions is directly related to the 
density, and gives similar results. It is clear from these flgures that there is a transition 
at the higher value of f3. We maintain that the lower-/3 line does not correspond to a 
thermodynamic phase transition. Indeed a similar line exists in the G model, shown in 
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Figure 6. Solutions of the renormalisation equation H17|) . The left-most line is 
consistent with the transition line conjectured in Machado et al ^J, whilst the right- 
most line is consistent with the transition line found using the condition Ai = A2, and 
shown in figure |31 The question remains as to which lines are true transition lines (see 
text). 

figure ini whilst it is known that there is no high-K transition hne. It is interesting to 
note that the hne does, however, seem to join up with the low-zt hne at the location 
of the tricritical 0-point. This is particularly interesting, since the equivalent line in 
the bond-interacting model also joins the line k = at a point consistent with the 

collapse transition (figure Ej), as noted by Machado et al [llj. It is less clear is whether 
this point is the same as the crystallisation point found by the condition Ai = A2 = 1, 
i.e. the termination point of the line k = K2{(3). We shall return to this issue later. 

In figure Uni the specific heat is plotted for /t = 0.5, calculated directly from the 
grand-canonical free energy using the fluctuation-dissipation theory: 

a = L/32(((e-M')-(e-/ip)'), (23) 

where e is the energy per site and /x is the chemical potential (/3/i = — log/t). 

The specific heat does not diverge, and has a broad peak, located at a lower value of 
j3 than the transition line for fixed k. This is reminiscent of the Berezinskii-Kosterlitz- 
Thouless (BKT) transition [2H], which is an infinite ordered phase transition, i.e. none 
of the derivatives of the free energy are singular. This identification is reinforced by 
the observation that the order parameter, 5p, is not related to a derivative of the free 
energy. It is known that walk models with geometrically frustrated interactions show 
BKT type transitions in the limit k oo|2n]- To the best of our knowledge, this would 
be the first time this has been observed at finite k. 

Figure ITD shows the phase diagram calculated using Nightingale's phenomenological 
renormalisation group method for odd lattice sizes. On the same diagram we have 




Figure 7. The order parameter, {pu — Pv\), suitable for the high-K transition (A) and 
its fluctuations (B) for k — 0.5. 



plotted the peaks of xsp^ which may be seen to converge nicely to the same line. These 
results agree with those found earlier for even lattice sizes, using the eigenvalue crossing 
method (see figure El). 

4. Locating the multicritical points 

In this section we shall try to pinpoint the location of the collapse transition and the 
crystallisation transition along the k = line in order to determine if these points 

are the same = P^), or if there are two distinct transitions ^ (32). 

The most common way of locating the collapse transition is to plot v estimates. For 
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Figure 8. The density, p, (A) and its fluctuations (B) for k = 0.5. 

T > Tcoii, the collapse-transition temperature, the estimates tend to the self-avoiding 
walk value, whilst for T < Tcou the estimates of u are expected to tend to 1/2; as 
n K,*, the correlation length may be identified with the linear size occupied by 
the walk. Since, at the first order transition, the walk fills the lattice to a finite 
density, its linear dimension scales as N^^'^, or \k — The value of u at the 

collapse transition is expected to take on an intermediate value. Crossings would be 
expected to converge to the point {Pcoii, ^coii) at the transition point. Additionally, ([T^ 
is simply a phenomenological renormalisation group equation, and so crossings of u may 
be identified with fixed points, corresponding to different universahty classes. 

Figure IT^ shows the estimates for u calculated using Nightingale phenomenological 
renormalisation, see (fT8|l . for the odd lattice sizes, for which Ai > A2 along the whole 
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Figure 9. Solutions of the rcnormalisation equation (|17|l for lattice sizes L = 3, L' = 5 
for the 0-point model, showing the existence of a high-density solution, which does not 
correspond to a known transition line. The dashed line corresponds to the transition 
line, and the solid line to the high density solution. The region where it may be 
expected to join the true transition line is consistent with the finite-size estimate for 
the 0-point. This seems to remain true as the lattice width is increased. 



0.5 - L=3 




0.5 1 1.5 

beta 



Figure 10. Plots for the specific heat for k = 0.5. The peaks of the specific heat 
converge to a finite value, at a value of f3 lower than the phase transition. The shape 
and behaviour are reminiscent of a KT type transition. 
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Figure 11. Proposed phase diagram calculated using the phenomenological 
renormalisation group equation 1)17(1 using using odd lattice widths. On the same 
diagram we plot the peaks of xsp- 

line. The estimates may clearly be seen to cross. There is also an intriguing series of 
shoulders, which seem to accumulate to form a fixed point with a value of u ~ 2/3, 
close to the branched polymer value for z/jHO]- This may be an indication that the walk 
forms branching structures before fully collapsing. This is more likely to be a crossover 
effect than a true phase transition, but this point warrants further investigation. 

The values calculated for the crossing points are given in tables ^ and |21 With so 
few data points it is not possible to extrapolate, particularly remembering that there 
are strong parity effects, however, looking at the odd lattice sizes, the transition is likely 
to be in the region of (3coi\ ~ 1-22 — 1.24. Values of u are consistent with a collapse 
transition in the 9-point universality class (2^3,5,7 = 0.525974, 1^5^,9 = 0.556653 and 
z/Q = 4/7 = 0.5714286). 

The location of the collapse transition can be estimated from the peaks of the 
response function for the corresponding order parameter. The density provides a good 
order parameter for the collapse transition, whilst 6p provides a good order parameter 
for the crystallisation transition. The plots of xp and xsp along the transition line are 
given in figure and the location of the corresponding peaks is given in tables 
and El The results derived from xsp indicate a transition which is certainly at a value of 
P > 1.27 whilst Xp indicates that the collapse transition occurs at a value of /3 < 1.25. 
This leads to the conclusion that there are at least two distinct multicritical points along 
the low-K transition line (k = 

We cross-check these results using additional estimates for the transition points. 
For the collapse transition we used phenomenological renormalisation using p (see (1201 - 
1221)); whilst for the crystallisation transition we looked for solutions to Ai = A2 = 1. 
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Figure 12. Estimates of v calculated using (|18|l . Using phenomenological 
renormalisation arguments, crossings may be identified with fixed points. Dashed lines 
indicate the three special values of v for the 8-point model: I'sAW — 3/4, vq — 4/7 
and J^densc = 0.5. We clearly see a fixed point corresponding to the pure self-avoiding 
walk model behaviour, as well as the correct dense walk limit as /3 is increased. Whilst 
the crossing between the L = 3,L' = 5 and the L = 5, L' = 7 line is far from 
h' — 4/7, the crossing between the L = 5, L' — 7 and L — 7, L' = 9 lines is rather 
close. Unfortunately, with so few estimates, it is not possible to extrapolate. Note the 
possible appearance of an additional solution with v « 2/3, see text for discussion. 



Table 1. Different estimates for the fugacity at the collapse transition, KcoU, using 
the different methods described in the text. 



L 


Nightingale 


f 


Xp 


Xsp 


Cv 


Ai = A2 = 1 


3 


0.268701 


0.258766 


0.37705 


0.30406 


0.31415 




4 


0.280168 


0.270168 


0.34815 


0.35443 


0.35781 


0.303687 


5 


0.277469 


0.267450 


0.31116 


0.26744 


0.28073 




6 






0.27017 


0.30231 


0.31440 


0.277016 


7 






0.28871 


0.25980 


0.27225 




8 






0.28516 


0.27886 


0.29573 


0.269245 


9 






0.28002 


0.25872 


0.26949 




oo 






0.271 ±0.003 


0.258 ±0.001 


0.267 ±0.001 





Plots of ^L,L+2 are shown in figure El and the estimates for the transition points using 
these two methods are, again, shown in tables Q and El The estimates for v found 
using the phenomenological renormalisation group are again consistent with a G-point 
transition. 
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Figure 13. Plots of (A) Xp, (B) Xip (C) the specific heat with n = 
the low-K-transition line. 
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Table 2. Different estimates for the value of P at the collapse transition, PcoU, using 
the different methods discussed in the text. 



L 


Nightingale 


f 


Xp 


X6p 


Cv 


Ai = A2 = 1 


3 


1.283504 


1.283169 


0.67946 


1.08935 


1.04256 


— 


4 


1.205824 


1.205822 


0.77646 


0.94286 


0.93014 


1.133416 


5 


1.240191 


1.240300 


0.98783 


1.24091 


1.17257 




6 






1.03850 


1.13837 


1.05522 


1.220231 


7 






1.11407 


1.27827 


1.21309 




8 






1.13936 


1.20315 


1.12213 


1.247185 


9 






1.16533 


1.28481 


1.22818 




oo 






1.23 ±0.02 


1.29 ±0.01 


1.24 ±0.01 
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Figure 14. ipL.L+2{P) plotted along the low-K-transition line. Again we see crossings 
corresponding to the self- avoiding walk fixed point, and the collapse transition. For the 
Q point we would expect a value of (/3 = —1/4. From the above figure the estimate of 
ip is closer to —0.11 (i^ w 0.53), although the value of ip does decrease as L is increased. 
Again with just two solutions, it is not possible to extrapolate. 



5. Discussion 



There are two basic phase diagrams proposed in the hterature for this model, both 
based on extended mean-field type calculations. By the mean-field nature of the 
calculations, they can only propose general features, and not details of critical behaviour. 
This is amplified in walk models, since the connectivities, and thus the single walk 
nature of the model, cannot be taken into account in the local description inherent 
in the methods used. However, our transfer matrix calculations tend to confirm the 
qualitative predictions made by Buzano and Pretti[lfl, rather than those made by [Hlini- 
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In summary, we propose the phase diagram shown quahtatively in figure |21 Moving 
along the self-avoiding walk line, the walk first collapses by the standard 0-point type 
transition, followed by a further transition where the corners are ejected from the walk, 
and the walk lines up with one of the two lattice directions. This latter transition is 
similar to the transition seen in the so called F-model of a gas of semi-rigid self-avoiding 
loops filling the lattice with a density equal to one (k oo). In the F-model the 
transition is in the BKT universality class. We have found evidence that the bond- 
interacting self-avoiding walk model also displays a BKT transition, but unusually this 
transition is extended into a line of transitions in the high-K phase. 

Whilst the 6-point model does not show a dense walk phase transition, other 
models do, notably the Hydrogen-bonding self-avoiding walk modelj^l, the Blote- 
Nienhuis 0{n) walk model [21] and the G-point model with the inclusion of a bending 
interaction jUlJ 1221 IHEl OH IHS]- These models all have one feature in common: they 
include interactions which are geometrically frustrated. It is expected that this 
frustration is seen when the density of the walk becomes finite, where the fractal (or 
Hausdorff) dimension is equal to the lattice dimension |22j. In the bond-interacting 
model the exponent found for the collapse transition is clearly larger than 1/2, hence 
the Hausdorff dimension is lower than the lattice dimension, and we could therefore 
expect the critical behaviour not to be modified by any lattice effects. In this case it is 
reasonable to expect the collapse transition to be in the universality class. 

It is interesting to note that in the three models in which the dense-phase transition 
line has been studied beyond mean- field, i.e. the Blote-Nienhuis 0{n) model, the 
Hydrogen-bonding model, and the bond-interacting model, all seem to have different 
classes of transition: for the Blote-Nienhuis 0{n) model the dense-phase transition line 
was in the Ising universality class |36j. in the Hydrogen model the transition was critical 
with a divergent specific heat with a value of u < and now the bond-interacting 
model with a BKT transition. At first sight it seems as if there are as many different 
thermodynamic behaviours as models one might define, but preliminary results for a 
Hydrogen-bonding self-avoiding walk, extended to include solvent quality effects (in the 
language of a polymer model) indicate that this model displays at least the latter two 
behaviours. 

Whilst transfer matrix calculations enable an exploration of the whole phase space, 
not accessible to Monte-Carlo simulation, the method is hampered by the small number 
of lattice widths available. Recently we introduced an extension to the CTMRG method 
for walk models 1^ . but unfortunately this method is not easy to implement for 
bond-interactions. Monte-Carlo simulation may be useful to study the limit ^ oo, 
enabling an independent study of the different transitions proposed. 
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